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ABSTRACT 

We generate new standard solar models using newly analyzed nuclear fusion cross sections and 
present results for helioseismic quantities and solar neutrino fluxes. We discuss the status of the 
solar abundance problem and investigate whether nonstandard solar models with accretion from the 
protoplanetary disk might alleviate the problem. We examine a broad range of possibilities, analyzing 
both metal-enriched and metal-depleted accretion models and exploring three scenarios for the timing 
of the accretion. Only partial solutions are found: one can bring either the depth of the convective zone 
or the surface helium abundance into agreement with helioseismic results, but not both simultaneously. 
In addition, detailed results for solar neutrino fluxes show that neutrinos are a competitive source 
of information about the solar core and can help constrain possible accretion histories of the Sun. 
Finally, we briefly discuss how measurements of solar neutrinos from the CN-cycle could shed light 
on the interaction between the early Sun and its protoplanetary disk. 

Subject headings: Sun: helioseismology - Sun: interior - Sun: abundances - neutrinos - accretion 



> 

O 



X 



1. INTRODUCTION 

Solar models have been very successful in describing 
the detailed structure of th e solar interior that emerged 
from helioseismic studies dChristensen-Dalsgaard et all 
[l^lBahcall et al.ll2001HTurck-Chieze et al.ll2001h . Fur- 
ther success came with the first direct measurement of 
the total 8 B neutri no flux by the SNO Collaboration 
(|Ahmad et al.l I2002D , which proved to be in excellent 
agreement with stan dard solar model (SSM) predictions 
(jBahcall et al.l [200Th . Solar models are sensitive to the 
details of their constitutive physics. Helioseismic predic- 
tions, for example, depend critically on radiative opaci- 
ties and microscopic diffusion, and are also sensitive to 
uncertainties in input parameters, such as abundances 
determined from absorption line analyses of the Sun's 
atmosphere. Predictions of solar neutrino fluxes are sen- 
sitive to the rates of key nuclear reactions. Solar models 
have played a mayor role not only in understanding the 
Sun, but also as a test bench for constituitive physics 
of stellar models in general. The further development 
of SSMs - both improved treatments of the physics and 
more accurate input parameters - as well as tests of the 
model against new solar observations are thus of funda- 
mental importance to stellar astrophysics. 

Over the last decade, the advent of 3D hy- 
drody namic models of near -s urface solar c onvec- 
tion (IStein fc Nordlundl Il99l lAsplund et al.l l2000l : 
iFrevtag et al.l2002HWedemever et al.l2004ft. a more care- 
ful se lection of spectral lines (e.e. lAUende Prieto et al.1 
12001 and, in some cases, relaxation of the assumption 
of local thermodynamic equilibrium in line formation 



(lAsplund et al.l l2004t iBergemann et al.1 12010D have re- 
sulted in a thorough reassessment of the photospheric 
solar a bundances. Th e most recent and complete revi- 
sion bv lAsplund et al.l (|2009t) . hereafter AGSS09, shows 
a reduction in the abundances of the volatile CNO ele- 
ments and Ne (—0.1-0.15 dex) w ith respect to older com- 
pilations of solar abund ances (|Grevesse fc Noels! Il99l 
iGrevesse &; Sauval 119981 hereafter GS98), as well as 
smaller changes for abundant refractory elements, e.g. Si 
and Fe. One simple figure-of-merit for comparisons is the 
photospheric total metal to hydrogen fraction, (Z/X) Q 
of -0.018 for AGSS09 versus -0.023 for GS98. The 
AGSS09 results, however, have not gone unchallenged. 
Based on 3D solar a tmosphere models c o mputed with 
the C Q5BOLD code (jFrevtag et al.l 120021) . ICaffau et al.l 
(|2010l and references there in) have derived best- value 
solar CNO abundances that are systematically higher 
than those in AGSS09 and associated uncertainties that 
are larger, so that the ranges are formally consistent with 
GS98 values. The differences between the two groups 
seem to arise not from the 3D hydrodynamical atmo- 
sphere models, but from the selection of spectral lines 
and assumptions in the line formation calculations. 

Nuclear reaction rates are another important input 
for SSMs. A very broad effort to critically evaluate pp 
chain and CN cycle rates wa s undertaken in 1 997-8 , 
culminating in Solar Fusion I lAdelberger et al.l (|1998l 
hereafter SFI). Since SFI a great deal of additional 
experimental a nd theoretical w o rk has been done, e.g. 
14 N(p , ^) 15 Q (iFormicola et all 12004 llmbriani et al l 
[2005t iRunkle et al.l 120051: IBemmerer et al.l I2006bt 
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iLemut et all l2006t IMarta et al.l 120081). 
(iSingh et aLll2004j: iBemmerer et all 



20071: IConfortola et al.1 "72007: 



3 He( 4 He,7) 7 Be 
2006at iGvurkv et all 
Brown et al.l 120071: 



di Leva et all 120091), a He( 3 He, 2p) i He (jBonetti et all 



1999ft. and 7 B e(p.7) 8 B (iHammache et al.l f l998. 



Strieder et ail 120011: Uunghans et all I2002T 



2001; 
2002 



Babv et all l2003al )bl: Uunghans et"alT 120101 ). The 



progress includes new facilities for measuring cross 
sections at low energies in or near the Gamow peak, 
concerted efforts by independent teams to address key 
uncertainties such as the S-factor for 7 Be(p, 7) 8 B, and 
the development of ah initio methods for theoretical 
extrapolations to solar energies. In order to assess 
the impact of this new work on the solar model and 
other applications, a second evaluation was recently 
undertaken, establishing new recommended values and 
uncertainties for pp ch ain and CN cycle cross sections 
(jAdelberger et al]|2010t hereafter Solar Fusion II; SFII). 
Here we present new solar models in which the SFII 
recommendations are adopted as new standard values, 
thereby providing the most current best values for SSM 
neutrino flux predictions and other solar properties. 

Standard solar models (SSMs) are calibrated against 
a present-day surface (Z/X) Q value: a standard as- 
sumption in SSMs is the equivalence between solar 
core primordial metal abundances and contemporary 
surface abundance^. Early work on SSMs incorpo- 
rating reduced CNO abundances (IBahcall et al.l l2004t 
IBasu fc Antiall2004l IMontalban et al.ll2004l) showed that 
such solar models are typically problematic: the loca- 
tion of the boundary of the convective envelope (Rcz), 
the sound speed and density profiles, and the present- 
day surface helium abundance (Is) predicted by SSMs 
using low (Z/X)q all differ from the values extracted 
from helios eismology by amoun ts beyond accepted un- 
certainties (Ba hcafi et ail 120061) . On the other hand, 
SSMs using higher (Z/X)q values, e.g. from GS98, show 
excellent agreement with all helioseismic constraints 
(iChristensen-palsgaard et al.lfl996t iSchlattl et all 119991 : 
IBahcall et al.ll200U l2005d) . The disagreement between 
'low-Z' solar model predictions and helioseismic con- 
straints, reminiscent of the now solved 'solar neutrino 
problem', is popularly known as the 'solar abundance 
problem'. 

Various attempts have been made to solve the solar 
abundance problem. For example, solar abundance de- 
term in ations have be e n ques ti oned (lAntia fe Basull2005l 
2006 1 : IBahcall et alj l2005bl; iDelahave fe PinsonneauTtl 
20061 iDelahave et al.ll2010l ) and the const itutive physics 
of solar models has been r evised (IBasu fc Antial 
2004 IMontalban et all 120041 : IBahcall et al.l 12005a: 
Guzik fc Mussackl 120101 ). If we assume the correct- 
ness of the low-Z metallicity determination for the Sun, 
there is not, at the moment, a satisfactory solution that 
brings the solar model into agreement with helioseis- 
mology. The only way to achieve agreement in a low- 
Z model comparable to that achieved with high metal- 
licity (GS98) models is through ad hoc adjustments, 
such as an increase in opacities to compensate for the 
lower metal abundance in the so lar radiative interior 
(|Christensen-Dalsgaard et al.ll2009l ). 

1 The equivalence is, of course, modulo the changes produced by 
microscopic diffusion. 



An al ternative pos si bility was considered by [Guzik 
(120061). iCastro et all (12007ft . and more recently by 
IGuzik fc Mussackl (|2010D~ These authors assumed that, 
soon after the Sun settled on the Main Sequence (MS), 
very metal-poor (or -free) matter was accreted onto the 
convective envelope, mixing with and diluting its primor- 
dial composition. Solar models were constructed with 
low-metallicity AGSS09-like convective envelopes, but 
with high-metallicity radiative interiors similar to GS98. 
It was hypothesized that such a model might have an in- 
ternal structure consistent with helioseismology, despite 
its low-metallicity photosphere. However, large discrep- 
ancies in model sound speed profiles and convective en- 
velope depths persisted. Additional modifications to the 
models (extended overshooting, turbulent mixing) were 
explored, but the difficulties remained, particularly the 
discrepant sound spee d profile below the convective en- 
velope (see Figure 1 in ICastro et al.lf2007h . 

While the possibility of solving the solar abundance 
problem through an inhomogeneous model is intriguing 
- we know metal inhomogeneities were established late 
in the evolution of the proto-planetary solar disk because 
a large excess of metal w as sequestered in the planets 
(jHaxton fc Sere nclli 2008) - it may be difficult to guess 
a priori the specific structure such a model should take. 
The work reported here is a first attempt to examine 
a broad class of possible nonhomogeneous solar mod- 
els that could be generated through accretion, which 
we model dynamically. We explore a range of accre- 
tion histories, varying the mass and composition of the 
gas as well as the epoch during which the accretion oc- 
curs. We study the detailed helioseismic properties of 
these models as well as their neutrino flux characteris- 
tics, to determine whether any of the models might be 
compatible with the various observational constraints. 
We also consider the implications of the models for fu- 
ture proposed neutrino tests of core composition, such 
as mig ht be possible through de tection of CN-cycle neu- 
trinos (|Haxton fc S erenclli 2008) - what classes of accre- 
tion models could be directly ruled out by such neutrino 
measurements? This paper establishes a baseline for fu- 
ture accretion studies that could be more directly tied 
to recent evidence for peculiar differences in the surface 
abundances of the Sun and similar stars (solar twins), 
and possible connections these differenc es may have to 
the c hemistry of proto-planetary disks (jMelendez et al.l 
2009). These additional studies will be presented else- 
where. 

The paper is laid out as follows. In § [2] we present 
updated calculations of SSMs for high and low metal- 
licity solar compositions, incorporating the most recent 
nuclear reaction rates. These models represent the base- 
lines against which models with accretion are later com- 
pared. In § [3] we describe our calculational procedures: 
the numerical algorithm and the parameter space the 
calculations span. Results for accreting solar models are 
presented in § [4] and discussed in § [5] Final remarks and 
conclusions are presented in § [6] 

2. STANDARD SOLAR MODELS 

The new SSM calculations presented here were done 
with GARSTEC (|Weiss fc Schlattlll200l) . including the 
up dates and mod i ficatio ns in the input physics described 
in iSerenelli et al.l ((2009). The new models, in addition, 
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adopt the nuclear reaction rates that were recently rec- 
ommended in SFII. 

SFII represents the first systematic evaluation of hy- 
droge n burning reacti ons since the SFI (lAdelberger et al.1 
119981 ) and NACRE (jAngulo et all 11999( 1 efforts. SFII 
provides a set of standard S-factors and uncertainties 
that reflect the progress made in laboratory and theo- 
retical nuclear astrophysics over the last decade. The 
team of evaluators carefully considered which data sets 
were sufficiently complete and well documented to be 
included; how to combine data sets that span different 
energy ranges and that might be somewhat inconsistent 
in their error-bar claims; and how to associate reason- 
able uncertainties to fitting functions (when the data sets 
span the region of astrophysical interest) or to theoret- 
ical extrapolations of data (when no data exist in the 
solar Gamow peak). The treatment of uncertainties was 
conservative, due to the increasing relevance of the SSM 
as a laboratory for weak interaction studies and stellar 
evolution theory. The principal changes between SFII 
and SFI are: 

o SFI identified 7 Bc(p,7) 8 B as the least well understood 
reaction of the pp chain, and called for additional mea- 
surements, particularly below 300 ke V. The SFI eval- 
uation was based on a single data set (jFilippone et al.l 

H983allb1 ) due to concern that other results were ei- 
ther too poorly documented or otherwise not suffi- 
ciently reliable for inclusion. The SFI recommended 
S-factor is Si 7 (0) = 19±| eV b. The value us ed in 
previous SSM calcul ations (ISerenelli et all 12009( 1 was 
S 17 (0) = 20.6 eV b (|Junghans et al J 120031) . The SFII 
eval uation benefited from four new, high-quality data 
sets (IHammache et a!lll99a[200TtlStrieder et allfe OOl: 

Uunghans et all I2002L 12003: iBabv et al l I2003allbl ) and 
yielded Si 7 (0)=20.8 ± 0.7 (expt) ± 1.4(theor). The 
uncertainty is dominated by the required theoretical 
extrapolation. The estimated total uncertainty at so- 
lar Gamow peak energies is 7.5%. 

o The S-factor for the controlling reaction of the CN 
cycle, 14 N(p,7) 15 0, was assigned a large and asym- 
metric uncertainty in SF I, Sn4(0) = 3.5 jl^ e 
b. A series of new measurements ranging from 
70 to 480 keV dFormicola et al.l 120041: llmbriani et all 
2001 IRunkle et all 120051: IBemmerer et al.l I2006bt 
Lemut et al.l 120061: flVIarta et all 120081 ) has led to a 



significant reduction in the best value and uncer- 
tainty: the SFII result, after summing all transitions, 
is Siu(0) = 1.66 ± 12 keV b. This val ue is close to 
Sn 4 (0) = 1.57 keV b lMarta et all (120081) adop ted in 
our previous SSM results ISerenelli et al.l ((20091) . The 
value and uncertainty are important to future plans 
to use CN neutrinos as a probe of the solar-core C 
and N abundances, as well as to age estimates for the 
galaxy's oldest stars. 

o Improved measurements were made of the reac- 
tions that govern the ppl/ppll+pplll branching ra- 
tio, 3 He( 3 He,pp) 4 He and 3 He( 4 He, 7 ) 7 Be. The SFI 
uncertainties for these reactions were 7.4% and 9.4%, 
respectively. At the time of SFI, the most recent 
measurement of the 3 He( 4 He,7) 7 Be cross section was 
ten years old. The new data employed in SFII 



for 3 Hc( 3 Hc,pp) 4 Hc include measurements that ex- 
tend to the lower edge of the solar Gamow peak 
(Bo netti et all 119991 ) . thus removing any dependence 
on theory apart from screeening corrections. The new 
data for 3 He( 4 H e,7) 7 Be inclu de a se ri es of high-quality 
measu r ements (|Singh et all 20041: IBemmerer et all 
2006 a! iGvurkv et al l 120071: IConfortola et al.l 120071 



Brown et al. 2007: Idi Leva et all 120091) at both high 



and low energies, spanning the range from 93 keV to 
2.51 MeV. The resulting SFII recommended S-factors 
are consistent with SFI values but significantly more 
precise, S 33 (0) = 5.21 ± 0.27 MeV b and S 34 (0) = 
0.56 ± 0.03 keV b. The central value of the SFII 
ratio S33(0)/S 3 4(0) is 8.7% lower than that of SFI, 
increasing somewhat ppll+pplll burning. Previously 
(ISerenelli et alll2009D. we had adopted S 34 (0) = 0.567 
keV b ICostantini et all ((20081 ) and S 33 (0) = 5.4 MeV 
b (SFI). 

o The S-factor for the driving reaction of the pp chain, 
p(p,e + ^ e )d, is taken from theory, as we have no di- 
rect means of measuring this reaction currently. The 
SFI value, Su(0) = 4.00 ■ 10~ 22 (l±0.007t^ 2 °) keV 
b, includes estimated uncertainties from input experi- 
mental quantities such as the axial- vector coupling qa 
and theoretical uncertainties reflecting possible vari- 
ations in potential-model calculations, including cor- 
rections from short-range physics and two-body cur- 
rents. The most significant development since SFI is 
the recognition that the overall nuclear physics un- 
certainty in Sn is governed by a single low-energy 
parameter which can be constrai ned by using the pre - 
cisely know tritium j3 decay rate ((Carlson et al.lfl99ll ). 
Consequently t hree distinct classes o f calculations, po- 
tential theory dSchiayilla et allll998D. hybrid effective 
field theory (EFT) dPark et all I2003D . and pionless 
EFT (iKong fe Ravendah]||200lHButler fe Chenl 120011 : 
lAndo et al.l 120081 ) . have now converged in their pre- 
dictions. The SFII recommended S-factor is Sn(0) = 
(4.01 ± 0.04) ■ 10 -22 keV b, and the uncertainty in 
the Gamow peak is estimated to be 0.9%. The cen- 
tral value i s alm ost 2% higher than what originally 
iPark et a l. (2003) recommended and was used in pre- 
vious SSM calculations. 

o The reaction 3 He(p,e + ;/ e ) 4 He is the source of the pp 
chain's most energetic neutrinos. The SFI S-factor 
best value is S hep (0) ~ 2.3- 10~ 20 keV b, but the evalu- 
ators concluded that no meaningful uncertainty could 
be assigned. They noted that the scatter of theoretical 
estimates corresponded to a factor of 2.5, up or down. 
After SFI th ere were significant effo rts in both poten- 
tial models dMarcucci et. al.1 120001) and hybrid EFT 
((Park et all 12003( 1 to improve the calculations and to 
assess possible sources of uncertainty. The two ap- 
proaches are in good agreement. The recommended 
SFII value is S hcp = (8.6 ± 2.6) • 10" 20 keV b. 

o New data f rom the LUNA coll aboration on 
15 N(p,7) 15 (jBemmerer et all 120091) tha t agree 
with an older data set from (Hebbardl (|1960l ) but not 
with iRolfs fe Rod"nevl (jl97J) led to the SFII con- 
clusion that some unidentified syste matic error may 
have affected the low-energy data of IRolfs fc Rodnevl 
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TABLE 1 

Main characteristics of 
standard solar models. 





GS98 


AGSS09 


(Z/X) s 


0.0229 


0.0178 


Z S 


0.0170 


0.0134 


Y S 


0.2429 


0.2319 




0.7124 


0.7231 


(Sc/c) 


0.0009 


0.0037 


z c 


0.0200 


0.0159 


Y c 


0.6333 


0.6222 


Zini 


0.0187 


0.0149 


Y ini 


0.2724 


0.2620 



(jl974l ). A significantly lower value for S7 15 (0) = 36 ±6 
keV b was recommended as an interim value, pending 
the outcome of very recent experiments now being 
analyzed. 

o The SFII ratio of pep to pp rates is 2.6% higher 
than the corresponding SFI ratio because the for- 
mer includes a calc ulation of radiative corrections 
(jKurvlov et alJl2003l ) while the later was evaluated at 
tree level. 

These new rates, and others from SFII not explicitly dis- 
cussed above, were incorporated in the new solar model 
calculations. 

All models in this work are required to reproduce the 
present-day (tq = 4.57 Gyr) solar luminosity, Lq = 
3.8418 x 10 33 erg s" 1 , and radius R Q = 6.9598 x 10 10 cm. 
In this section we present two SSMs that differ in the 
present-day surface metal-to-hydrogen abundance ratios. 
One model is based on the GS98 solar composition, 
(Z/X) s = 0.0229, the other on AGSS09, (Z/X)* = 
0.0178; in both cases the adopted scale is meteoritiqj. 

In Table [1] we list the most relevant parameters for 
these SSMs. The models are labeled according to the 
adopted solar compositions. The results differ little from 
those of earlier mod els with the same resp ective compo- 
sitions, presented in lSerenelli et al. (2009, Table 2). The 
SFII reaction rates do not introduce noticeable changes 
in the initial composition of the models, denoted here by 
the initial metal, Z lnl , and helium, i^ n ;, mass fractions. 
Consequently, the surface and central mass fractions of 
metals, Z§ and Zq, and of heliu m, Y"s and Yp, a re als o 
almost unchanged from those of iSerenelli et al.l {2009). 
The depth of the convective envelope Rcz is slightly 
deeper in the models with the SFII rates and the av- 
erage fractional rms deviation between the solar and the 
model sound speed profiles, (6c/ c), is also slightly bet- 
ter in the new models. The changes are, however, small 
compared to model uncertainties, so the new models are 
effectively fully consistent with previous ones. 

Given the similaritie s betw een the new models and 
those of ISerenelli et al.1 ()2009l ). the reader is referred to 
that source for a more detailed discussion. Here, we fo- 
cus on manifestations of the solar abundance problem, 
the primary motivation for the studies presented here. 
There are three key indicators of the abundance prob- 
lem, the depth of the convective envelope, the surface 

Except for volatile elements; see lAsplund ef al. (2009) for de- 
tails. 



TABLE 2 

Solar neutrino fluxes. 



v flux GS98 AGSS09 Solar a 



PP 


5.98(1 ±0.006) 


6.03(1 ± 0.006) 


6.05(lt° ;2g) 


pep 


1.44(1 ±0.012) 


1.47(1 ± 0.012) 




hep 


8.04(1 ± 0.30) 


8.31(1 ±0.30) 


i8(i±8:ib 


7 Be 


5.00(1 ± 0.07) 


4.56(1 ±0.07) 


4-82(iig;8S) 


8 B 


5.58(1 ± 0.14) 


4.59(1 ±0.14) 


5.00(1 ±0.03) 


13 N 


2.96(1 ± 0.14) 


2.17(1 ±0.14) 


< 6.7 


15Q 


2.23(1 ± 0.15) 


1.56(1 ±0.15) 


< 3.2 


17 F 


5.52(1 ± 0.17) 


3.40(1 ±0.16) 


< 59. 


x 2 /pa gr 


3.5/90% 


3.4/90% 





Note. - Neutrino fluxes are given in units of 

10 10 (pp), 10 9 ( 7 Be), 10 8 (pep, 13 N, ls O), 10 6 ( 8 B, 17 F) and 
10 3 (hep) cm -2 s . Asymmetric uncertainties have been av- 
eraged. 

a Solar neutrino fluxes in ferred from all available neutrino data 
(Borcxino collaboration 2011). 

helium abundance, and the sound speed profile. The 
AGSS09 model produces a convective zone that is too 
shallow and a surface helium abundance that is too low 
(see Table The entries in the table should be com- 
pared to the h elioseismic values, Rcz = 0.713 ± 0.001 
ilBasu fc Antidll997l ) and 0.2485 ±0.0035 (|Basu fc Antial 
2004) , respectively. If helioseis mic uncertainties ar e com- 
bined with SSM uncertainties (jBahcall et al. 2006, equa- 
tions 32 and 35), then AGSS09 model predictions differ 
from solar values at 3a for both Rcz and is- ln con- 
trast, GS98 model predictions are in excellent agreement 
with the parameters extracted from helioseismology. Al- 
though more difficult to quantify, the difference between 
the AGSS09 and solar sound speed profiles, determined 
from (5c/ c), is about four times larger than that ob- 
tained with the higher metallicity GS98 model. These 
specific discrepancies for the AGSS09 abundance set are 
the quantitative manifestations of the solar abundance 
problem mentioned in the introduction. 

The new SFII nuclear reaction rates do alter the pre- 
dicted solar neutrino fluxes, particularly the fluxes asso- 
ciated with the ppIII chain and CN-cycle, mechanisms 
for solar hydrogen burning that are relatively unimpor- 
tant energetically. W ith respect to the rates used in 
ISerenelli et al.l (|2009f ) . the most important SFII changes 
are the 2% increase in Sn(0), the 3% decrease in 833(0), 
and the 5% increase in Si 14(0). Other changes affecting 
predicted neutrino fluxes include the reduced proton cap- 
ture cross section on 15 N and the increase in the pep/ pp 
rate ratio. Model predictions for the neutrino fluxes and 
associated uncertainties are presented in TableOfor both 
the AGSS09 and GS98 compositions. The most signifi- 
cant changes are a 5% decrease in the predicted 8 B flux 
primarily because of the increase in Sn and the increase 
in the 13 N flux due to the larger central abundance of C. 
The increase in C is a consequence of the lower SFII value 
for 15 N(p,7) 16 0, a reaction that competes with the CN 
I cycle reaction 15 N(p,a) 12 C and allows mass to flow out 
of the CN I cycle into CN II. We discuss the parametric 
dependence of the fluxes on nuclear cross sections and 
other SSM input parameters elsewhere, including some 
of the resulting neutrino tests that might be made of 
solar composition and S-factors (|Serenelli et al.ll201lT ) . 
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Table [5] also includes the updated solar neutrino fluxes 
inferred from all available neutrino data. The analysis in- 
cludes the recent more precise, 7 Be measurement, which 
is the main change wi t h resp e ct to previous analysis 
(IBahcall fc Pena-Garavl [2001 IPena-Garav fc Serenellil 
2008; Gonzal ez-Garcfa et al.ll201Qfl/ Details of the anal- 
ysis will be shown elsewhere (Borexin o collaboration! 
120111 ). In order to compare the SSM predictions with the 
fluxes inferred from neutrino data, we use th e \ 2 func- 
tion defined in iGonzalez-Garci'a et al.l ((2010), with up- 
dated errors and correlations. We find Xqs98 = ^-5 an d 
XAGSS09 = 3 - 4 > leading in both cases to PgIqs.agssoq = 
90%. The new fusion cross sections from SFII and the 
new Borexino results lead to both models predicting solar 
neutrino fluxes in excellent agreement with inferred ones; 
from Table [2] it can be seen that solar 7 Be and 8 B fluxes 
are intermediate between SSM predictions for GS98 and 
AGSS09 compositions. Currently, solar neutrinos can 
not differentiate between different solar compositions. 

3. SOLAR MODELS WITH ACCRETION: 
COMPUTATIONAL DETAILS 

3.1. The algorithm 

Briefly, the construction of an SSM requires finding 
two parameters that determine the initial composition 
of the model, Yini and Z- ln i (Xi n j = 1 — Yj n i — Zim), and 
the mixing length parameter «mlt- The three parame- 
ters are determined iteratively by requiring the model to 
reproduce the present-day solar luminosity, radius and 
surface metal-to-hydrogen abundance ratio. We remind 
the reader that the choice of a solar composition deter- 
mines not only the present-day (Z/X) s value, but also 
the relative abundances of metals within the bulk metal- 
licity. In what follows, all models are computed using 
the AGSS09 solar composition. 

The algorithm used to construct SSMs can be easily 
modified to compute nonstandard models with accre- 
tion, once the properties of the accreted matter (mass 
and composition) and the timescales and accretion rates 
are specified. Below, and in the remainder of the pa- 
per, we use the following nomenclature: M ac is the mass 
of the accreted material, Z ac , A ac , and Y ac are, respec- 
tively, its metal, hydrogen and helium mass fractions, 
T ac ,i and r aCjC are the initial and ending times of the ac- 
cretion process. Schematically, the algorithm consists of 
the following steps: 

1. An initial pre- main sequence (PMS) model of mass 
M = Mq — M ac is constructed using initial guesses 
for Y^, Z- mi , and «mlt- 

2. The model is evolved at constant mass until r = 
r aCi i. The accretion then begins at a constant rate 
M, continuing until r = r ac c , when M = Mq. 
Z ac is a parameter of the model, while A ac and 
Y ac are determined by the conditions A ac /Y ac = 
Ai n i/Yni and 1 = A ac + Y ac + Z ac . That is, we 
assume that the hydrogen-to-helium ratio of the 
accreted material is identical to that of the initial 
solar material. 

3. From r = T aCi0 onward the model is evolved at 
constant mass until r = r Q . Model properties are 



checked against present-day solar values, and the 
differences are used to improve the initial guesses 
for Yini, Z; n i, and «mlt- The process is repeated 
until the solar surface properties are reproduced 
to the desired accuracy (typically 1 part in 10 4 or 
better). 

In the above scheme the free parameters of the model 
are the same as those of SSMs, namely Y] n i, Z iln , and 
c^mlt- But they now depend on properties of the ac- 
creted material, particularly its composition, because the 
mixture of initial and accreted matter that forms the 
present-day solar surface must satisfy the (Z/X) s con- 
straint. 

So far, we have assumed that the parameters charac- 
terizing the accretion process are known. We now discuss 
how we selected the ranges for those parameters. 

3.2. Models 

The properties of the resulting solar models will clearly 
depend on the parameters that describe both the proper- 
ties of the accreted mater and the timing and duration of 
the accretion process. Our goal was to cover a relatively 
large parameter space, in order to better define the range 
of possible solar model outcomes. Below we describe our 
choices. 

3.2.1. Metallicity of the accreted mass 

Intuitively, a hypothetical solution to the solar abun- 
dance problem would seem to require that metal-poor or 
even metal-free material is accreted onto the Sun. This 
has b een the assumption in the models explored previ- 
ously (|Castro et al.ll2057t iGuzik fc Mussack1[20l0l ). Here 
we explore metal-free, metal-poor, and metal-rich accre- 
tion. Surprisingly, as described in § [4j we find that metal- 
rich accretion can lead to partial improvements in the he- 
lioseismic properties of some models. We explore Z ac val- 
ues ranging from to 0.03, the latter being roughly twice 
the solar surface metallicity. The relative abundances 
of metals is not varied, but is fixed to AGSS09 values. 
This is an assumption that may be a reasonable start- 
ing point for initial explorations, but probably should be 
relaxed in future calculations given that any mechanism 
for segregating metals would likely reflect condensation 
temperatures and other element-specific chemistry. 

3.2.2. Accreted mass 

We consider accretion of mass from the protoplanetary 
disk with a composition different from that of the unpro- 
cessed material of the proto-solar nebula. The amount 
of such material the Sun might have accreted in its early 
evolutionary phases in not well constrained, as estimates 
of both the mass of the proto-planetary disk and the frac- 
tion of that mass that remained in the solar system vary. 
Assuming that initially the proto-solar cloud was well 
mixed, one can set a lower limit to the mass of the pro- 
toplanetary disk, the so-called minimum-mass solar neb- 
ula (MMSN) limit, from estimates of the metal content 
of the planets. The calculation, howev er, still depends 
on model assumptio ns. Early estimates (jWeidcnsc hillingl 
[l977t IHavashi|[l98ll) give M M msn ~ 0.013 - 0.07 M Q , 
assuming planets were formed at their present locations. 
More sophisticated models that tak e into accoun t mi- 
gration give somewhat larger values. iDeschl (|2007ft finds 
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Mmmsn ~ 0.17 M Q but with a disk that is being strongly 
photo-evaporated, i.e., most of its mass migrates out- 
wards, not towards the Sun. 

Based on these considerations, we explore values of 
ilf ac hi the range from 3 • 10 -5 to 0.06 M Q , or equiv- 
alently, from 10 to 2 ■ 1O 4 M0. We stress that this is 
the total accreted mass, not the metal mass. The up- 
per end of the range thus comes within a factor of 2-3 of 
most estimates of Mmmsn- Given that models of mas- 
sive disks wi th steep density profiles lead to substantial 
(|Deschll2007[ ) evaporation, 0.06 M Q is probably a safe up- 
per bound on the amount of mass that could have been 
chemically processed in the protoplanetary disk and later 
accreted onto the Sun. In terms of the total accreted 
metal, our choices for M ac and Z ac yield a range of to 
600 M®. 

3.2.3. Accretion timescales and mass accretion rate 




Age [Myr] 

Fig. 1. — Time dependence of Mqe during early evolutionary 
phases of a solar model; the shaded area depicts the convective 
zone. M/Mq = 1 defines the surface of the model. 

The time when accretion takes place is of particular im- 
portance. This is clear from inspection of Figure[TJ where 
the evolution of the mass of the convective envelope, 
Mce, during early SSM evolution is shown. Initially, 
the model is fully convective on the Hayashi track. As it 
evolves towards the Zero Age Main Sequence (ZAMS), a 
radiative core develops and convection quickly retreats. 
During these early phases, the e-folding time for the evo- 
lution of Mce is about 8 Myr. At about 30 Myr, the star 
reaches the ZAMS and further evolutionary changes, in- 
cluding the extension of the convective envelope, occur 
on very long (nuclear) timescales. For a fixed amount 
of accreted material, the mass of the convective envelope 
at the time of accretion determines the degree of dilution 
(or enrichment). These considerations lead us to identify 
three paradigmatic accretion scenarios: 

o early accretion: Accretion starts when the convective 
envelope still encompasses most of the star; accreted 
material is then highly diluted. We take T aCj i = 5 Myr 
as a representative time for the early accretion sce- 
nario. In typical models, the convective envelope con- 
tains 70 to 80% of the total solar mass at this time. 

o intermediate accretion: Accretion starts when the con- 
vective envelope has receded towards the surface and 



contains about 20 to 30% of the stellar mass. Ac- 
creted material is strongly diluted in the envelope but 
the effects of accretion in the final composition of the 
model are clearly noticeable. We take r aCi i = 15 Myr 
as representative of this case; 

late accretion: Accretion starts when the model is very 
close to the ZAMS, i.e. when timescales for struc- 
tural changes begin to be dominated by nuclear burn- 
ing timescales. The convective envelope is mature, 
close to its final MS size, and changes very little after- 
wards. Its mass is comparable to some of the largest 
M ac values we consider. Accordingly, the final struc- 
ture of models with late accretion is strongly affected 
by the properties of the accreted material. We take 
Tac,i = 30 Myr as representative of late accretion. 

Accretion could occur at even later times. However, 
given that Mce changes very little after the first 30 Myr 
of evolution, such models should not differ significantly 
from the late accretion scenario described above. In 
terms of involved timescales, the la te accretion scenario 
is e quivalent to that investiga ted by iCastro et al.l ()2007f ) 
and iGuzik fe Mussackl (|2010f ). who assumed accretion 
took place once the Sun had settled on the MS. 

Computationally, the accretion scenarios are defined in 
terms of the time at which accretion takes place. Phys- 
ically, however, the relevant characteristic is the thick- 
ness of the convective envelope. Therefore, as long as 
accretion occurs sufficiently early that element diffusion 
effects are still negligible, the thickness of the convec- 
tive envelope is the only physically relevant parame- 
ter. For this reason, even if the young Sun avoided 
a ful l y convective Hayashi ph ase (|Wuchterl fc Klessenl 
l200lt iBaraffe fc Chabrierl l2010h and consequently devel- 
oped a thin convective e nvelope after only a few Myr 
(|Baraffe fc Chabrierl[20ll . our intermediate and late ac- 
cretion results remain valid, constraining these alterna- 
tive evolutionary tracks. 

The duration of the accretion phase (AT ac ) is a second 
important parameter, particularly for the early and inter- 
mediate cases, where Mce evolves rapidly. Our standard 
choice is Ar ac = 10 Myr. There are two reasons for this 
choice. One is simplicity: this choice corresponds to qui- 
escent accretion rates well below 10 -7 M Q yr _1 , allowing 
us to neglect the energy deposition on the outer layers 
of the star and associated effe cts that could alter final 
internal structure of the model (jBaraffe et aL 2009). We 
assume that the central star is almost completely assem- 
bled when accretion from the protoplanetary disk, that is 
accretion of chemically processed material, begins. The 
second reason is the correspondence to typical estimates 
of the planetary formation time, usually in the range of 
a few to 10 Myr. However, to test the impact of differ- 
ent accretion durations, we have computed models with 
more rapid accretion, Ar ac = 1 Myr. 

In all cases considered, accretion is assumed to occur 
at a constant mass rate M = M ac /Ar ac . 

3.3. Grid of models 

The ranges considered for M ac and Z ac were explored 
by stepping through a 10 x 10 grid, with quadratic spac- 
ings. Therefore, for each value of r aCj i and Ar ac , 100 
models were generated. This basic grid is depicted in 
Figure [H 
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Fig. 2. — The mass and metallicity grid for the model calcula- 
tions. 

4. SOLAR MODELS WITH ACCRETION: RESULTS 

Here we present results for the three accretion scenar- 
ios discussed above, early, intermediate, and late. We 
first focus on the relationship between the altered solar 
composition and helioseismic observables, then later dis- 
cuss the effects on neutrino fluxes. 

4.1. Early accretion 
4.1.1. Composition and helioseismic quantities 

Model compositions and resulting helioseismic quan- 
tities are shown in Figure for five selected model se- 
quences of constant Z ac , from metal- free to super-solar 
metallicity. Panels a and b give the initial composition 
of the models, Y n i and Z lv ^ respectively, i.e. the compo- 
sition before accretion takes place. Panel c shows the 
present-day surface metallicity Zs- The lower panels 
show selected solar properties that can be directly tested 
with helioseismology, Rcz in panel d, {5c/ c) in panel e, 
and Ys m panel /. 

We first examine the behavior of the initial compo- 
sition of the models. Since all models have to satisfy 
(Z/X)s = 0.0178, metal-poor accretion models lead to 
higher Z- ln i values. A higher metallicity core implies an 
increase in the radiative opacity and molecular weight, 
i.e., larger central temperatures. However, models are 
constrained by L Q as well. So, an increase in nuclear 
burning caused by higher central temperatures has to be 
compensated by limiting the amount of fuel available, 
that is increasing Yni- In panel a we include the pri- 
mordial helium mass fraction from Standard Big Bang 
Nucleosynthesis (SBBN), Y P = 0.2486 ± 0.0006, using 
th e latest WMAP res ults for the baryon-to-photon rato 
■q (jLarson et al.ll2010j ). 

The effects of accretion on metallicity profiles are 
shown in Figure @] for three models with M ac = 0.06. 
The models are identified in panel e by the open circles. 
For comparison, we included the results for the SSMs 
presented in § [2] The strong gradient in the metallicity 
profile between ~ 0.17 and 0.35 R© seen for the metal- 
rich and metal-free models is the result of mixing be- 
tween the accreted and the primordial material, while 
the convective zone is receding towards the surface. As 
discussed above, the surface metallicity is very similar in 
all models and nearly equal to that of the AGSS09 SSM. 



Model values for Rcz, {5c/ c), and Yg, quantities tightly 
constrained by helioseismology, are shown in panels d, e, 
and / of Figure 03 respectively. Panels d and / also in- 
clude the helioseismic values for Rcz and Ys, plotted as 
horizontal solid lines, and their la uncertainties, shown 
as dotted lines. We first discuss the behavior of Rcz as 
a function of Z ac . The outer layers (above 0.8 M Q , see 
Figure [4j of models evolved with metal-poor accretion 
have lower metallicities. Additionally, lower Zg requires 
lower Xg- B oth effects tend to l ower the radiative opac- 
ity at Rqz (Bahc all et al J 120041 ). producing a shallower 
convective envelope. The opposite is true for metal-rich 
accretion. Therefore, the behavior of Rcz seems to argue 
against metal- free (or metal-poor) accretion. In contrast, 
metal-rich accretion leads to deeper convective envelopes, 
closer to the helioseismic value. The trends in the behav- 
ior of the sound speed profiles are quite similar, as one 
sees from panel e. There the GS98 SSM result is shown 
as a horizontal line. The behavior of {5c/ c) and Rcz are 
very well correlated because the bulk of the discrepancy 
in the sound speed profile comes from the mismatch be- 
tween the solar and model Rcz values. To show this, we 
plot in Figure [5] the relative sound speed difference as 
a function of radius for the same models shown in Fig- 
ure 2] and identified in panel e with open circles. The 
bump below the convective zone decreases for models 
with Rcz close to the solar value (vertical dotted line). 
This occurs for metal-rich accretion models. The SSM 
results for AGSS09 and GS98 compositions are included 
for comparison. 

Unfortunately, the behavior of Ys is qualitatively op- 
posite that of Rcz and {5c/ c). Metal-poor accretion 
leads to higher Y~s va lues as a result of larger Yni 
(jSerenelli fc Basul fcOlOh improving the agreement with 
the helioseismic determination. In contrast, metal-rich 
accretion leads to lower, more discrepant, Ys values. 
Therefore, helium abundance considerations favor metal- 
poor accretion. 

These results suggest that it may not be possible to 
construct solar models with early accretion in which all 
predicted helioseismic variables agree with solar values 
extracted from helioseismic observations. While agree- 
ment for Rcz and {5c/c) can be improved with metal- 
rich accretion, the accompanying lower Ys values increase 
discrepancies for that variable. Similarly, metal-poor ac- 
cretion helps with Ys, but increases the discrepancies in 
Rcz and {5c/ c). 

4.1.2. Neutrino fluxes 

Solar neutrino fluxes place a second set of con- 
straints on solar interior properties. We will com- 
pare model predictions to the flux values determined 
from a global analysis of so lar neutrino experiments 
(|Borexino collaboration! 1201 ID that are given in the 
last column of Table [2j The authors considered all 
available data from solar neutrino experiments, in- 
cluding new determi nations of the 8 B flux by SNO 
(|Aharmim et al.ll2010D and Borexino ([Bellini et al.ll201Cft 
and the most recent measurement of the 7 Be flux by 
Borex ino (jArpesella et al.l 120081 : iBorexino collaboration! 

1ml). 

The general dependence of solar model fluxes on ac- 
cretion simply reflects the initial metallicity, which de- 
termines the temperature in the Sun's energy-producing 
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Fig. 3. — Composition and helioseismic quantities of sequences of accreting solar models with T aCj j = 5 Myr and Ar ac = 10 Myr (early 
accretion). Results are shown for five different values of Z ac , as labeled in panel a. 
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Fig. 4. — Final metallicity profile as a function of radius, for three 
accreting models with the Z ac values indicated in the plot and for 
M ac = 0.06 Mq. For comparison, the SSM metallicity profiles are 
also shown (see §0- 

core. Fluxes that depend more strongly on temperature 
( 8 B, 13 N, 15 0, 17 F and, to a lesser extent, 7 Be) increase 
with increasing core metallicity (metal-poor accretion). 
However, because the models are constrained by L = Lq, 
neutrino fluxes less depende nt on core temperature (pp, 
pep and hep) must decrease (|Bahcall fc Ulmerlll996l ). In 
our accretion models the dependence of each neutrino 
flux on solar core metallicity can be characterized by the 
partial derivatives 



dlog$(z^ 

8 log Zini 



(1) 



We determine Oi(z/j), for each flux i, by performing a joint 
linear fit to all early accretion models. The fluxes most 
strongly correlated with core temperature yield a( 7 Be) = 
0.34, a( 8 B) = 0.82, a( 13 N+ 15 0) = 1.58. For pp and pep 
fluxes we find a(pp) = —0.036, a(pep) = —0.16. In what 
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Fig. 5. — Sound speed profiles for the models of Figure [4] The 
vertical dotted line marks the boundary of the convective envelope. 



follows, the sum of the fluxes 13 N + 15 O will be denoted 
CN. 

Results for the various neutrino fluxes are shown in 
Figure [6] The metallicity sequences are the same as in 
Figure [3] Panels a and b show the 8 B and 7 Be fluxes, 
which have been measured directly in experiments (the 
7 Be flux can be extracted from Borexino results given 
other measurements of neutrino oscillation parameters). 
The solar fluxes (third column, Table [2]) are shown with 
horizontal black lines, including la uncertainties. Models 
with metal-rich accretion tend to deviate from the solar 
values, particularly for the precisely known 8 B flux due to 
its sensitivity to core temperature. As mentioned before, 
the greatest sensitivity to accretion is exhibited by the 
combined CN flux (panel c) because, in addition to its de- 
pendence on core temperature, this flux depends linearly 
on the total C+N abundance in the solar core. Unfortu- 
nately, the current determination of the combined flux is 
only an upper limit, $(CN) < 9.9 x 10 8 cm _2 s _1 , and 
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Fig. 6. — Neutrino fluxes for the accreting solar models of Figure [51 
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Fig. 7. — Composition and helioseismic quantities of sequences of accreting solar models with r ac i = 15 Myr and Ar ac = 10 Myr 
(intermediate accretion). Results are shown for five different values of Z ac , as labeled in panel a. 



does not allow to put any constraints on the solar core 
metallicity. 

Panels d and e show the strongly correlated pp and 
pep fluxes. The luminosity constraint and the 7 Be mea- 
surement by Borexino combine to produce the small un- 
certainties in these fluxes (about 1% and 1.5%, respec- 
tively). As with $( 8 B), particularly the pep flux some- 
what disfavors metal-rich accretion models, especially 
those with large M ac values. 

As an indicator of the overall agreement between model 
and solar fluxes, we calculated, for each model, the 
X 2 function (see § The results are shown in panel 



/. The horizontal dashed-dotted line shows the 68.3% 
probability band. As expected from results for individ- 
ual fluxes, metal-poor accretion models lead to a small 
overall improvement in the agreement with experimen- 
tal fluxes, mostly driven by improvements in the agree- 
ment with $( 8 B) and $(pep). Note that the sequence 
with Z ac = 0.0130 shows almost no changes in the pre- 
dicted neutrino fluxes and \ 2 values: this value of Z ac 
is nearly equal to the AGSS09 Zi n i value, so that the 
accreted material and initial solar material have nearly 
identical compositions. We caution that these neutrino 
flux changes should not be overinterpretted, as we have 
seen that models with metal-poor accretion tend to do 
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poorly on Rqz and (5c/ c) and therefore their internal 
structure shows differences with respect to the Sun. 

4.2. Intermediate accretion 
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Fig. 8. — Metallicity profiles for selected models (open circles in 
panel e in Figure 0. SSM results are also shown, as in Figure [4] 



4.2.1. Composition and helioseismic quantities 

The results for Yi nl , Z- ln i, and Zg are shown in panels 
a, b, and c, respectively, of Figure [3 The qualitative 
trends track those we observed in early accretion mod- 
els: metal-poor accretion leads to higher Z ini and YJni 
and, because of the lower X; n ;, lower Z$. Metal-rich ac- 
cretion produces the opposite behavior. Metal-rich ac- 
cretion for Z ac = 0.024 and 0.030 and with accreted 
masses above 0.05 and 0.027 M , respectively, lead to 
Y[ n i lower than the SBBN primordial value Yp . The same 
models lead to Z- m i values that are low when compared 
to m etallicities observed in the solar neighborhood (see 
e.g. iPrzvbilla et al.l [20081 ). Compared to early accre- 
tion models, the consequences of accretion are amplified 
in the intermediate models because the convective en- 
velopes are shallower at the time of accretion. The large 
contrast between Zi n i and Zg for intermediate accretion 
models is apparent from Figure [SI where the metallic- 
ity profiles for the three models identified with open cir- 
cles in Figure [7] are shown. The partially mixed region 
extends roughly between 0.35 and 0.55 R©, depending 
on the model. Profiles of the GS98 and AGSS09 SSMs 
are also shown. Metal-rich accretion models, as stated 
above, lead to under-metallic cores. The consequences 
for resulting model sound-speed profiles are significant, 
as discussed below. 

Results for the helioseismic quantities Rcz, (5c/ 'c), and 
Y$ are shown in panels d, e, and / of Figure respec- 
tively. The qualitative behavior of Rqz and Yg is quite 
similar to that found in the early accretion models. In- 
terestingly, models with metal-rich accretion can lead to 
Rcz m excellent agreement with the helioseismic value. 
However, the same models yield Yg values that are low 
compared to helioseismic values, as well as Y[ n [ values be- 
low the SBBN value. For Z ac = 0.030, the highest M ac 
values considered, the convective envelopes are deeper 
than allowed by helioseismology. This is reflected in the 
behavior of (5c/ c). For the Z ac = 0.03 sequence, (5c/ c) 
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Fig. 9. — Sound speed profiles for the models of Figure l8l 

has a minimum around M ac = 0.027 M©. The agreement 
is better than with the AGSS09 SSM, but still a factor 
of two worse than with the GS98 SSM. The agreement 
deteriorates for larger M ac values, as the convective en- 
velope becomes too deep. Figure Ogives the sound speed 
profiles for the same models that were used in Figure [8j 
The sound speed profile for the model with metal-free 
accretion (black solid line) is significantly degraded com- 
pared to that of AGSS09 SSM, because the convective 
envelope is too shallow; the bump below the convective 
zone is very large and extends quite deep, to ~ 0.35 R©. 
The difficulty exists despite the similarity between the 
metallicity of this model and that of AGSS09 SSM over 
a substantial region, 0.45-0.7 R©, below their convective 
envelopes (Figure O for M > 0.9M Q ). For Z ac = 0.024 
and M ac = 0.06M© (green line-dotted line), the model 
predicts the correct i?cz, so that the bump below the 
convective envelope is almost gone. At small r, the sound 
speed difference oscillates around zero, and the overall 
(5c/ c) is twice as large as for the GS98 SSM (shown as the 
grey dashed line). Finally, the model with Z ac = 0.030, 
which we noted had too deep a convective zone, exhibits 
a large bump in the sound speed profile of opposite sign 
below the convective zone, as well as significant discrep- 
ancies at smaller radii due to its under-metallic interior. 

4.2.2. Neutrino fluxes 

Relative to early accretion models, the neutrino fluxes 
from intermediate accretion models exhibit slightly more 
sensitivity to Z ini : a( 7 Be) = 0.52, a( 8 B) = 1.12, 
a(CN) = 1.77, a(pp) = -0.045, and a(pep) = -0.14. 
Flux results are shown in Figure [TOl The higher sensitiv- 
ity of temperature-dependent fluxes, particularly $( 8 B), 
$( 7 Be), and the combined $(CN), reflects the decreased 
dilution of accreted material. Panel a shows that $( 8 B), 
in models with metal-rich accretion, differs markedly 
from the solar value (horizontal black line). Conversely, 
metal-poor accretion modestly improves the agreement 
with experiment for almost all the range of M ac we con- 
sidered. Metal-free accretion can also lead to an im- 
proved $( 8 B), but only for smaller M ac values. Oth- 
erwise, for models with M ac > 0.04M Q , the predicted 
$( 8 B) is too large. Because 7 Be (panel 6) is less sensitive 
to core temperature changes (and less precisely known), 
significant discrepancies arise only for metal-rich accre- 
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tion at large M ac values. Panel c shows the sensitivity of 
the combined CN flux to the accretion, but this flux is 
too poorly constrained experimentally to allow any con- 
clusions to be drawn. 

The results for the pp and pep fluxes (panels d and e) , 
though dependent on the luminosity constraint, strongly 
suggest that metal-rich accretion can be excluded for 
Af ac > 0.03 M©. However, metal- free and metal-poor 
accretion, over the range of M ac values we considered, 
lead to models that are perfectly consistent with the pp 
and pep solar values. This conclusion is consistent with 
the one we drew from examining the 7 Be flux. 

The global \ 2 m panel / shows that, as in the early 
accretion case, the metal-free and metal-poor accretion 
scenarios can slightly improve the overall agreement be- 
tween the models and the the solar neutrino fluxes. 
It should be recalled once more, however, that both 
AGSS09 and GS98 SSMs give excellent agreement with 
experimentally determined sloar neutrino fluxes (Ta- 
ble The horizontal lines in the plot are the 68.4%, 
5%, and 0.1% probability levels. Comparing the x 2 val- 
ues of the different Z ac sequences with results for the in- 
dividual neutrino fluxes, one finds that the 8 B flux, with 
its small uncertainty and strong temperature sensitiv- 
ity, dominates the \ 2 ■ F° r example, although the model 
with Z ac = 0.00 and M ac = 0.06M Q yields nearly perfect 
agreement with the Be, pp, and pep fluxes, a relatively 
high x 2 is obtained because the model overestimates the 
8 B neutrino flux. 

4.3. Late accretion 
4.3.1. Composition and helioseismic quantities 

In a late accretion model the accretion occurs when the 
Sun is close to the MS with a well-developed radiative 
core and consequently a thin convective envelope. The 
effects of the accretion on the convective envelope are 
therefore maximal. 

The composition and helioseismic properties of the 



models are summarized in Figure [TT1 The sequences for 
Z ac = 0.024 and 0.030 terminate at M ac = 0.02 and 
0.012M Q respectively. As is evident from the plot of Z; n ; 
in panel 6, larger values of M ac are incompatible with 
modern Sun's (Z/X)g: the needed initial metallicity is 
driven to negative values. Even for smaller, allowed M ac 
values, metal-rich accretion leads to very low Yi n j values, 
well below the SBBN value, as shown in panel a. 

The helioseismic quantities i?cz, (6c/ c), and Is are 
plotted in panels d, e, and /, respectively. The metal-free 
and metal-poor accretion models (solid black and blue 
short dashed lines) yield shallow convective envelopes. 
Also, for relatively low M ac values (M ac < 0.015 M Q ), 
the overall sound speed profiles, as estimated from (5c/ c) 
and Ys, are somewhat improved. These models have ini- 
tial compositions and thus core compositions very simi- 
lar to that of the GS98 SSM. However, as the models do 
poorly on Rqz , discrepancies in the sound speed profiles 
remain. Models with larger M ac values at low Z ac fail to 
reproduce the observed Rcz and (6c/ c), and also over- 
estimate Ys- We discuss metal-poor accretion further in 
§ 15. 3[ where we compare our results with previous work. 

In contrast to early and intermediate accretion mod- 
els, sound-speed profiles for late, metal-rich models do 
not improve, regardless of the M ac value. The accre- 
tion's dramatic effects on surface composition lead to 
unacceptably low values Z- 1TI \ values even for small M ac 
values. The resulting low-metallicity cores are not com- 
patible with the radiative-zone sound speed profiles. Yg 
is also too low for these models. 

4.3.2. Neutrino fluxes 

Neutrino fluxes are quite sensitive to parameter varia- 
tions in late accretion models, responding to the large 
changes in Z- lri \. The power-law exponents we derive 
are: a( 7 Be) = 0.66, a( 8 B) = 1.43, a(CN) = 2.08, 
a(pp) = -0.081, and a(pep) = -0.18. 

The individual fluxes are given in Figure 1121 The 
agreement between model and observed fluxes deterio- 
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Fig. 12. — Neutrino fluxes for the accreting solar models of Figure [TT] 



rates in both metal-free and metal-rich accretion mod- 
els. This is mostly driven by the 8 B flux (panel a), but 
the other fluxes with relatively small uncertainties ( 7 Be, 
pp, and pep) also tend to trend away from observation. 
Thus the neutrino fluxes limit the accreted composition 
to near-SSM values, independent of the M ac value. The 
slightly metal-poor Z ac = 0.0093 sequence is an inter- 
esting one, yielding for M ac ~ 0.01M Q a somewhat im- 
proved fit to the observed fluxes. In this model (5c/ c) 
and Ys are also somewhat improved, while the connec- 
tive envelope is only slightly shallower than that of the 
AGSS09 SSM. 



5. DISCUSSION 

We have investigated whether the solar abundance 
problem can be solved or alleviated by constructing so- 
lar models with different interior and envelope composi- 
tions. Rather than postulating an ad hoc modern Sun 
with an inhomogeneous structure, we instead generate 
that structure through a physical mechanism, accretion 
of metal-poor or metal-rich material from the protoplan- 
etary disk. This mechanism is well motivated both be- 
cause significant accretion is seen in T Tauri stars, which 
may resemble the early Sun, and because we see clear 
evidence in our solar system of the the large-scale segre- 
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gation of metal that accompanied planetary formation. 
The chemical processes that concentrated dust and ice at 
the disk's midplane, and the physical process of planetary 
formation that scoured those metals from the disk, lock- 
ing them into the planets, could have combined to create 
conditions for either metal-poor creation accretion (depo- 
sition of metal-depleted gas from the disk's outer surface 
layers) or metal-rich accretion (proto-planets spiraling 
into the Sun) after the Sun had developed its radiative 
core. In this way the composition of the Sun's convec- 
tive zone could have been altered, becoming chemically 
distinct from the Sun's radiative core, which would have 
a composition quite close to that of the solar system's 
primordial gas. 

By exploring accretion as the candidate mechanism 
responsible for an inhomogeneous modern Sun, we can 
generate nonstandard models as a function of physical 
parameters, such as the mass and composition of the ac- 
creted gas and the time of accretion. The accretion pa- 
rameters, and thus the resulting models, are constrained 
to reproduce the composition of the modern solar sur- 
face, taken to be AGSS09 in our models. We can then 
test these models by comparing their properties against 
contemporary solar observations done through either he- 
lioseismology or solar neutrino spectroscopy. The goal is 
to determine whether one can account for modern obser- 
vations that appear to be contradictory in the context 
of the SSM - the solar abundance problem - through 
a plausible physical mechanism that could have negated 
one of the SSM's key assumptions, a homogeneous ZAMS 
Sun. 

The possibility that metal-poor accretion could 
be responsible for the so lar abundanc e prob l em has 
motiv at ed earlier work (iCastro et al.l 120071: iGuzik I 
20061: lHaxton fc Serenellil 120081: iNordlundl 120091 : 
Guzik fc Mussackl l2010h . Our work extends these 
earlier efforts by exploring a broad range of accretion 
histories, masses, and compositions, and by testing 
the resulting models for consistency with modern 
observations of the Sun's surface and interior. However, 
there remain unexplored degrees of freedom: we have 
assumed a constant mass transfer rate, a fixed duration 
of the accretion, and fixed composition of the accreted 
material. We discuss our findings and their implications 
below. 

5.1. Helioseismic quantities 

The most apparent manifestations of the solar abun- 
dance problem are the inconsistencies between obser- 
vations and SSM solar surface helium abundances (too 
low), convective envelope depths (too shallow), and the 
solar sound speed profiles, particularly in the upper ra- 
diative zone. The main helioseismic results for our up- 
dated SSMs can be found in Table [1] 

How do the model predictions change if we postulate 
metal-rich or metal-poor accretion? Let us first consider 
results for Rcz and Is for the early, intermediate, and 
late accretion scenarios, shown in panels d and / of Fig- 
ures [H and [HJ respectively. Qualitatively, metal-poor 
accretion leads to higher Yg values and shallower con- 
vective envelopes with respect to the SSM. The opposite 
happens for metal-rich accretion. Therefore, accretion 
can be invoked to construct models that match the he- 
lioseismic values of either Yg (metal-poor) or Rcz (metal- 
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Fig. 13. — Top panel: average mis sound speed for the extended 
sequence of early accretion models, for Z ac = 0.030. Best model, 
for M ac = 0.090 Mq, is marked with an open circle. Bottom 
panel: the red long-dashed line is the sound speed profile for the 
M ac = 0.090 M model. The GS98 (grey solid line) and AGSS09 
(grey short-dashed) SSM profiles are included for comparison. 



rich), but at the expense of enlarging the discrepancy in 
the other quantity. No model simultaneously improves 
both Ys and Rcz- 

The discrepancies we observe between solar and model 
sound speed profiles are correlated with the discrepan- 
cies in Rcz- Better agreement in Rcz leads in general 
to smaller (Sc/c) values. This is a consequence of {8c/ c) 
being dominated by the large bump in the sound speed 
difference right below the convective envelope that is ac- 
centuated whenever there is a significant mismatch be- 
tween the solar and model values for Rcz- The best 
results for intermediate (late) accretion lead to {8c/ c) 
being about two (three) times the value for the GS98 
SSM. Interestingly, the best results for intermediate ac- 
cretion are from a metal-rich case (panel e in Figures [7]) 
and for late accretion are from a metal-poor case (panel 
e in Figures II 1 [) . The latter provides an exception to 
the general rule of correlated {8c/ c) and i?cz, stemming 
from the sensitivity of Z- m { in late accretion models to 
the properties of the accreted material. We discuss this 
further in § 15.31 

For the early accretion scenario, the results displayed 
in panel e of Figure [3] for metal-rich accretion suggest 
that further improvement in {Sc/c) might be achieved 
if M ac > 0.06 Mq. To test this, we computed mod- 
els with Z ac = 0.030 and Af ac = 0.075, 0.090, 0.110, 
and 0.125 Mq. The results for {8c/c) are given in 
Figure QS] (top panel). A minimum is found around 
— 0.090 M Q , yielding a value for {8c/ c) close to 
that of the GS98 SSM (solid black horizontal line). For 
accretion masses greater than M ac — 0.090 M©, {8c/ c) 
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begins to rise again. In the bottom panel the sound speed 
profile for the M ac — 0.090 M Q model is shown (red long- 
dashed line). It is similar to that of the GS98 SSM (grey 
solid line), except close to the Sun's center, where the 
metallicity and helium abundances in the model are too 
low. The AGSS09 SSM sound speed profile is also shown. 
For this model Rqz = 0.716 R©, in reasonable agreement 
with the helioseismic value 0.713 ±0.001 R©. However, 
following the general relation between Rcz and Y"s previ- 
ously described, this model has Is = 0.225, further from 
the helioseismic value, Y$ = 0.2485 ± 0.0035, than the 
AGSS09 SSM value (Y s = 0.2319, the limit of no accre- 
tion). Therefore, larger accretion masses do not offer a 
global solution to the solar abundance problem. 

All calculations presented in Section [4] have been com- 
puted assuming a fixed duration for the accretion phase 
of Ar ac = 10 Myr. Estimated life times of protoplan- 
etary disks are about a few Myr (Haisch et al.1 120011 : 
iWilliams fc Ciezal 1201 ID so that longer timescales for 
accretion seem unlikely. On the other hand, shorter 
timescales are possible, perhaps even favored. We there- 
fore repeated the calculations from § |4] using Ar ac = 
1 Myr, to test the consequences of shorter accretion 
times. As one might expect intuitively from an exam- 
ination of Figure [U the effect of a reduced Ar ac for given 
Z&c, M ac , and r aC! ; is analogous to enhancing the di- 
lution factor of the accreted material: that is, for two 
otherwise equivalent accretion scenarios with the same 
starting times, the one with the shorter duration will de- 
posit material into a larger mean convection-zone mass 
during the accretion, and thus that material will ex- 
perience more dilution. For example, results for (T ac ,i, 
Ar ac ) = (15 Myr,l Myr) are intermediate between those 
for (5 Myr, 10 Myr ) and (15 Myr, 10 Myr), shown in Fig- 
ures 03 and respectively. For late accretion, e.g., after 
the model has settled on the MS, the duration of the ac- 
cretion phase is irrelevant, provided it is short compared 
to gravitational settling timescales: at late times, there 
is little evolution in either the depth or the mass of the 
convective zone. Thus we conclude that changes in the 
duration of the accretion will not alter our basic conclu- 
sion: no combination of the accretion variables M ac , Z ac 
and r aC! i will lead to a model in which all helioseismic 
predictions are improved. 

5.2. Solar neutrinos 

We have compared model solar neutrino fluxes with 
solar fluxes (column 3 in Table [5]). The latter have been 
derived from a combined analysis of all solar neutrino 
experiments with the addition of the so lar luminosity 
constraintQ (B orexino co llaborat ionll2011[ ). The analysis 
includes the CN-cycle fluxes, although only upper lim- 
its are currently available. We find, using the new nu- 
clear reaction rates from SFII and the newest Borexin o 
results for the 7 Be flux (|Borexino c ollabora tion! [2011D . 
that the GS98 and AGSS09 SSMs are both in excellent 
agreement with solar neutrino data, producing compara- 
ble fits. These results are summarized in the last row of 
Table H 

Among the four well-determined fluxes, the 8 B neu- 
trinos are the most sensitive to accretion, responding to 

3 The luminosity constraint equates the solar luminosity to the 
rate of energy generation in the core, assuming a steady-state Sun. 



the variations in initial core metallicity that occur be- 
cause the accretion model is solved under the constraint 
of a fixed final convective-zone metallicity-to-hydrogen 
ratio. The 8 B solar flux is also the most certain, now 
determined to ~ 3%. Moreover, $( 8 B) provides infor- 
mation complementary to helioseismology because the 
production region, R < 0.1 R©, is essentially inaccessible 
to p-mode study. In fact, we find that all models with 
metal-rich accretion, because of their lower Z- lnl values, 
produce too few 8 B neutrinos. Results are shown in panel 
a of Figures [6j QUI and H2] This suggests that metal-rich 
accretion would create a conflict smaller than but rem- 
iniscent of the old solar neutrino problem, particularly 
for the intermediate and late accretion scenarios. In con- 
trast, models with metal-poor accretion, for an appropri- 
ate choice of M ac that depends on the accretion scenario 
(early, intermediate, or late), can be brought into excel- 
lent agreement with the experimental fluxes. The high 
sensitivity of the 8 B flux to metallicity leads to M ac val- 
ues as low as 0.01 M Q being disfavored for the extreme 
case of late metal-free accretion. 

The 7 Be, pp, and pep fluxes, are now also determined 
to a high level of precision, but are less sensitive to the 
core conditions (see Figures [H QUI and [T2|) . With respect 
to the AGSS09 SSM, a moderate amount of metal-poor 
accretion can lead to improvements in both the 8 B and 
7 Be fluxes, as discussed below. A slightly more metal- 
rich core seems preferred by these fluxes. 

We have done a global analysis of neutrino fluxes, with 
the results of the \ 2 calculations given in panel / of Fig- 
ures HI [TUl and Q21 In the case of early accretion, only 
the most metal-rich models with massive accretion show 
significant disagreement with solar data. The neutrino 
flux responses to accretion are more sensitive in the in- 
termediate and late scenarios, however. Figure [14] gives 
the resulting x 2 surfaces in the M ac — Z ac plane for these 
scenarios. In each panel the best-fit model is denoted by 
a red diamond. Solid lines depict the 68.3%, 10%, and 
1% probability contours, and dotted lines the contours of 
probability for the AGSS09 and GS98 SSMs (90%) (Ta- 
ble [2]). In the late accretion case, the top- right corner 
is empty because no solar models satisfying our require- 
ments can be constructed for such M ac — Z ac combina- 
tions. These plots define rather clearly the portions of the 
M ac — Z&c parameter space that are consistent with solar 
neutrino data. For intermediate accretion, all M ac val- 
ues considered for Z ac < 0.015 yield excellent fits, com- 
parable to SSM results. Improvements are possible but 
minor because, as we have noted before, AGSS09 SSM 
neutrino flux predictions are in very nice agreement with 
current solar neutrino data. We conclude that solar neu- 
trino data are primarily of value in ruling out accretion 
scenarios that give unacceptable fluxes (as opposed to 
providing positive evidence for accretion scenarios that 
marginally improve neutrino flux fits). As the accretion 
becomes more metal-rich, a larger range of high-M ac ac- 
cretion scenarios is excluded. In the case of late accre- 
tion, large fractions of the M ac -Z ac plane are significantly 
disfavored by solar neutrino data at both the metal-rich 
and metal-poor extremes. The full M ac range is allowed 
only for accretion models with Z ac w 0.015: these are 
cases close to the SSM. 

We close this section on solar neutrinos with a short 
discussion of the potential role of neutrino fluxes from the 
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Fig. 14. — Results for the combined x 2 global solar neutrino analysis for the intermediate and late accretion scenarios. Equal probability 
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the AGSS09 and GS98 SSM fits (Table [2} are overlaid, to indicate regions where accretion models do better or worse than the SSMs. 



CN cycle. These fluxes are very sensitive to metallicity 
because they respond in two ways: like the 8 B neutri- 
nos, the CN neutrinos respond sharply to environmental 
effects associated with metallicity - changes in the radia- 
tive opacities and the mean molecular weight that alter 
the temperature of the solar core. But the CN fluxes have 
an additional linear dependence on the core abundance 
of the CN elements, as C and N cat alyze the hydrogen- 
burnin g reactions of the CN I cycle. lHaxton fc Sereneliil 
(2008) have shown that a future precise measurement 
of $( 13 N) and $( 15 0) could be used to determine the 
solar core C+N abundance to ~ 10%, including all so- 
lar model and neutr ino parameter uncertainties (see also 
iSerenelli et al.ll2~011l ) . A direct determination of the cen- 
tral C+N abundance would be very valuable, as a com- 
parison to surface abundances would constitute a direct 
test of the SSM hypothesis of a homogeneous ZAMS Sun. 
In our models we find that the central-to-surface metal- 
licity ratio, normalized to the SSM value to account for 
the effects of microscopic diffusion, deviates from unity 
by more than 10% in many of our models, particularly 
for the intermediate and late accretion scenarios, as illus- 
trated in Figure [15l Thus a high-quality CN flux mea- 
surement could help rule out a large portion of the po- 
tential parameter space of accretion models, making the 
measurement an important diagnostics of the early Sun 
and its proto-planetary disk. Of course, this is in the 
context of the assumptions we have made in our accre- 
tion scenarios, which include a uniform rate of accretion 
of material with a metallicity related to AGSS09 by a 
simple scale factor, onto a proto-Sun with a metallicity 
that is also scaled to AGSS09. As disk condensation tem- 
peratures vary from element to element, there is ample 
motivation for relaxing these assumptions. (More general 
scenarios in which such assumptions are reconsidered will 
be discussed elsewhere.) Despite this caveat, we believe 
that the 13 N and 15 O neutrino fluxes are potentially very 
useful as probes of solar core metallicity and, therefore, 



potentially of the chemistry of the proto-solar disk. This 
chemistry sets the conditions for the formation of the so- 
lar system's diverse planets, and, as we have speculated, 
may have also altered the composition of the solar sur- 
face. Given that there exist very few quantitative probes 
of the processes that governed solar system formation, it 
is important to pursue any experiment that can directly 
constrain phenomena like accretion during the formation 
period. 

5.3. Comparison with previous works 

IGuzikl 12001: IGuzik fc Mussackl (|20loh and 
iCastro et al.l (|20071 ) have studied accretion of metal-poor 
material as a possible solution to the solar abundance 
problem. The authors constructed solar models in 
which the interior metallicity is high (around the GS98 
value ) and the enve l ope m eta llicity i s low, consistent 
with lAsplund etldl (|2005l) . IGuzik I (|2006t) assumes 
that a mass of 0.02 M Q is accreted during a period 
of 36 Myr, onto a ZAMS stellar model of 0.98 M . 
ICastro et al.l (j2007l ) takes a simplified approach in which 
the metallicity of the convective envelope is simply 
scaled down by a given factor (0.5) when the solar 
model is 74 Myr old. In both treatments the accretion is 
assumed to occur after the contraction phase to the MS 
is over, when the convective envelope is thin and very 
close to its present-day extension. In this aspect, their 
models resemble our late accretion scenario. 

In terms of Z^, two of our la te-accre t ion m odels 
closely resemble th e cases studied IGuzik I ([2006) and 
ICastro et all (|200l . These are our models with Z ac — 
0.0093 and with M ac = 0.013 and 0.018 M Q , respec- 
tively (see Figures Ql] and [12]). The sound speed pro- 
files of these two models, shown in Figure HH match 
very well that of the GS98 SSM below R/Rq = 0.5, 
but the bump below the convective zone remains quite 
prominent. Therefore {8c/ c) shows only modest improve- 
ment over the AGSS09 SSM value. These results corre- 
spond to Figure 12 of IGuzik I (g006) and to Fi gure 1 of 
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Fig. 15. — Ratio of the central to surface metallicity in selected 
models. Panels a, b, and c refer to early, intermediate and late 
accretion scenarios respectively. Metallicity sequences are shown 
for Z ac = 0, 0.0059, 0.0130, 0.0240, 0.0300 (top to bottom in each 
panel). 
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Fig. 16. — Sound speed profiles for the models with M ac = 0.013 
(black solid line) and 0.018 Mq (red dashed-dotted line), Z ac = 
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ICastro et all (|2007f) . Our models also yield R C z = 0.725 
and 0.727 R0, respectively, so the agreement with he- 
lioseismology in this quantity is slightly worse than for 
the AGSS09 SSM (0.723 R w ). These ob s ervat ions are in 
qualitative agreement with lCastro et all (|2007l ). who also 
find that the convective zone gets shallower (0.732 in- 
stead of 0.730 R Q ) in the model with accretion. iGuzik I 



(2006) , on the other hand, finds improvement in Rcz 
with respect to the model with no accretion. The rea- 
son for the qualitatively different behavior is not clear, 
but plausibly could be connected to the assumed ef- 
fects of the accretion on the convective zone's hydro- 
gen and helium abundances. Finally, our models give 
Ys = 0.244 and 0.250, i n good agreem ent with the he- 
lioseismic value, whereas iGuzik I (|2006T ) and ICastro et al.l 

(2007) find Y s = 0.240, marginally compatible with the 
seismic result. The small numerical discrepancies may 
be due to the use of somewhat different solar composi- 
tions in the three sets of calculations. Yet, overall the 
resulting models are in reasonable agreement, given the 
differences in the manner accretion was implemented in 
the three studies. 

6. SUMMARY AND CONCLUSIONS 

In this study we have explored the possibility that an 
accretion phase, occurring after the Sun developed a ra- 
diative core, could resolve, or at least mitigate, the solar 
abundance problem. The motivation for considering this 
possibility is the observation that large-scale metal seg- 
regation occurred in the proto-planetary disk, connected 
with the formation of the planets. Extending previous 
studies, we considered both metal-poor and metal-rich 
accretion and also allowed accretion to occur at different 
evolutionary stages of the young Sun. We also considered 
the constraints on models from both helioseismology and 
neutrino measurements, while earlier studies focused on 
the former. 

We began by updating the SSM, adopting the nuclear 
reaction rates from the very recent evaluation of SFII. 
The new SSMs for ASSS09 and GS98 abundances then 
were used as benchmark calculations against which the 
effects of accretion could be measured. The new mod- 
els differ only modes tly from others recently calculated 
(|Serenelli et al.l f2009h . The new solar models with the 
lower metallicity AGSS09 composition remain in con- 
flict with helioseismology. The newly recommended nu- 
clear fusion rates do impact predicted SSM solar neu- 
trino fluxes. The higher metallicity GS98 SSM remains 
in excellent agreement with solar neutrino data, but 
now the lower metallicity AGSS09 SSM shows an equiv- 
alent level of agreement. This also partly the result 
of the most recent 7 Be flux measurement by Borexino 
(jBorexino collaborationll2011[ ) . Thus solar neutrino spec- 
troscopy (and the associated laboratory astrophysics) has 
not yet reach the level of precision necessary to distin- 
guish between the AGSS09 and GS98 SSMs. 

The main conclusion we draw from our study of accre- 
tion is that, for the class of models we explored, no sce- 
nario was found to provide a satisfactory solution of the 
solar abundance problem. Here "satisfactory solution" 
is defined as a model with AGSS09 surface abundances 
that agrees with the known helioseismic properties as 
well as the GS98 SSM. We find that metal-rich accre- 
tion can bring the depth of the convective envelope into 
agreement with the seismic value, but that the resulting 
surface helium abundances in such models are always too 
low - in fact, in some of these models the initial helium 
abundance falls below the primordial SBBN value. Con- 
versely, some metal-poor accretion models bring surface 
helium abundance into a very nice agreement with the 
solar value, but at the cost of a convective envelope that 
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is even shallower than that of the AGSS09 SSM. 

Problems in the sound speed profile are generally 
strongly correlated with the problems in the depth of 
the convective envelope. The better (worse) the latter, 
the better (worse) the former. An exception to this gen- 
eral rule is found for late accretion models in which about 
0.015 Mq of metal-poor material is accreted. In this case, 
though the model's convective envelope is too shallow, 
there is some improvement in the sound speed profile. 
However, a significant bump in the sound speed differ- 
ence (model vs. heliosesimology) remains in the region 
right below the convective envelope. These models have 
surface helium abundances and neutrino fluxes that agree 
with solar values. 

Current solar neutrino data constrain solar core prop- 
erties. Our global comparison of model and solar fluxes 
rules out both metal-poor and metal-rich accretion for 
M ac > 0.01 M in the late accretion scenario. 

The results of our accretion study provide further mo- 
tivation for the proposal to exploit future measurements 
of CN neut rinos as a test of the C+ N abundance of the 
solar core (jHaxton fe Serenellil 12008). We find that the 
ratio between the surface and central solar metallicity in 
accretion models deviates from unity (gravitational set- 
tling taken into account) by more than 10% for most 
accretion scenarios. Such a measurement of the solar 
core bulk abundance of C+N does not directly constrain 
other metals and thus the overall metallicity: despite the 
simplifying assumptions made in our accretion model, 
there is no a priori reason to believe that the composi- 
tion of accreted gas is related to the primordial gas by 
a simple scale factor. Still, any quantitative abundance 
determination in the solar interior would be a major step 
forward. 



In the past, accretion of metal- rich material was dis- 
cusse d as a possible soluti on to the solar neutrino prob- 
lem ([Winnick et al.l I2002D . More recently, a new low- 
metallicity solar composition has generated concern that 
the Sun's surface might have been al t ered in an episode o f 
metal-poor accretion (jGuzik 1 120061 : iCastro et all 12007). 
Additional motivations for continuing the study of accre- 
tion onto the early Sun can be found in suggestions that 
the mechanism generating metal-poor or metal-rich ac- 
cretion might be connected to the dynamics of the proto- 
planetary disk. The possibility of such a connection has 
been raised in several papers: in the observation that 
solar abundance problem involves a deficit of convective- 
zone metal co mparable to the excess of m etal sequestered 
in the planets ( Hax ton fe Serenellill2008D: in recent mod- 
cls of the assembly of low mass stars (|Baraffe fc Chabrierl 
l2010f ) ; and in claims that the Sun has an abundance pat- 
tern of metals that correlates with metal condensation 
temperatures, suggesting a connection with disk chem- 
istry and the formation of rocky planets ([Mclcndc z et al.l 
2009). Work on accretion scenarios motivated by these 
observations is ongoing and will be reported elsewhere. 
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